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Abstract 

Emergent pattern formation in self-propelled particle (SPP) systems is extensively studied be- 
cause it addresses a range of swarming phenomena which occur without leadership. Here we present 
a dynamic SPP model in which a sensory blind zone is introduced into each particle's zone of in- 
teraction. Using numerical simulations we discovered that the degradation of milling patterns with 
increasing blind zone ranges undergoes two distinct transitions, including a new, spatially non- 
homogeneous transition that involves cessation of particles' motion caused by broken symmetries 
in their interaction fields. Our results also show the necessity of nearly complete panoramic sensory 
ability for milling behavior to emerge in dynamic SPP models, suggesting a possible relationship 
between collective behavior and sensory systems of biological organisms. 
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Self-organization and pattern formation in self-propelled particle (SPP) systems has been 
a topic of great interest in theoretical physics, mathematical biology and computational 
science It is well understood that the emergence of cohesive swarming motions requires 

neither leaders nor globally enforced organizational principles. Various SPP models have 
been used to explore stability and phase transitions of swarming patterns in response to 
varying noise levels ^,3,0] and other control parameters 3, Q,[8, 9|, as well as to characterize 



distinct regimes in the parameter space 



fly, 



10 



, [ll|. Effort has also been allotted for 



addressing biological questions concerning swarming behavior 



14, Q, 3 and 



for d esig ning non-trivial swarming patterns from combinations of different kinetic parameter 



sets 



17|. 



SPP models may be classified into two distinct categories: kinematic and dynamic 
Kinematic SPP models typically assume that each particle maintains a speed and orientation 
in accord with its local neighbors |l| . In these models, particles cannot halt because they are 
supplied with a minimal or constant absolute velocity y, Q, Q • Kinematic models have 
been used for computational modeling of collective behavior of constantly moving groups, 
such as bird flocks, often implementing empirically constructed complex, spatially discrete 
interaction zones and behavioral rules that reflect perceptional or locomotive properties of 
the species being modeled Q, Q, 1^, Q. On the other hand, dynamic SPP models 
describe the motion of particles using differential equations based on Newtonian mechanics 
that involve self-propulsion and pairwise attraction/repulsion forces jj]. It is known that 
such models may robustly form coherent millin g p atterns from initially random conditions 



even without explicit alignment rules 



y, y, 



10 



11| . Dynamic models have been used 



for both analytical and numerical research on collective behavior of interacting particles in 
general, with minimal complexity assumed in particles' intrinsic behaviors. 

Here we consider a new dynamic SPP model in which a sensory blind zone is introduced 
into each particle's zone of interaction. Although the assumption of sensory blind zones has 



been widely adopted in kinematic SPP models ^2 



13 



14 



impti' 

HE 



Il6l | , it has not been considered 



within a dynamic framework. We specifically examine the effect of the sensory blind zones 
on coherent milling behavior in dynamic SPP models. In doing so, we discovered a novel 
transition that occurred with an increasing range of blind zones, and found the system to 
be highly sensitive to this type of perturbation. 

Our model describes the movement, within an open, two-dimensional, continuous space. 



2 



of self-propelled particles driven by soft-core interactions whose dynamics are given by 

= {a-(3\v,\')v,-VU,{x,), (2) 
Ui{x) = J2u{\x-Xj\), (3) 

«(r) = ae-"/''- - Cae-"/'» (r > 0), (4) 

where Xi and Vi are the position and the velocity of the i-th particle {i = 1 . . . N), respec- 
tively; m the unit mass of one particle; a and (3 the coefficients of propulsion and friction, 
respectively; Ui{x) the interaction potential surface for the i-th particle; u{r) the pairwise 
interaction potential function (Fig. [1] (a)); Cr and Ca the amplitudes of repulsive and at- 
tractive pairwise interaction potentials, respectively; and Z,. and la the characteristic ranges 
of repulsive and attractive pairwise interaction potentials, respectively. Eq. ([2]) includes a 
velocity-dependent locomotory term and an interaction term achieved through a generalized 
Morse pairwise interaction potential. For a,f3 > 0, particles will rapidly approach equilib- 



rium velocity of magnitude v^q = yct/(3 and the system will converge toward a structure 
for which total dissipation is zero and particles are driven only by conserved forces . This 
rule set has been employed, with some mathematical variation, by many previous studies 



In this study, the shape of the pairwise interaction potential falls wit 



hin the 

biologically relevant regime defined as Cr/Ca > 1 and Ir/la < 1, as described by In 
the biologically relevant regime, individuals tend to move toward other individuals that are 
further from, and away from individuals that are closer than, some critical distance from 
themselves. This rule is generally applicable to the kinetics of many different biological 



species and natural systems 
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We compare two experimental parameters in the above model: the magnitude of stochas- 
tic force (noise) 7 and the range of sensory blind zones 6, the former analyzed in 5| and 
the latter our original extension. Sensory blind zones are incorporated into the design of 
each particle to mimic the abilities of anisotropic sensory systems observed in nature, such 
as vision. A sensory blind zone is assumed to exist for each particle with an angular range 
^ in a direction opposite to the direction of forward motion (Fig. [1] (b)). 

The inclusion of these parameters introduces discrete events into the model, i.e., abrupt 
changes of velocity by stochastic force and entry and exit of other particles into/out of 
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5,,,= {2,3,4,5,7} 

FIG. 1: Model assumptions used in the SPP model employed by this study, (a) Shape of the 
pairwise interaction potential function u{r) defined by Eq. ([4]), where r is the distance between 
two particles. This study explores the biologically relevant regime of parameter settings, in which 
particles will accelerate away from neighbors who are closer than, and toward neighbors further 
than, the equilibrium distance r^q = ^"-^^ log 0^ (~ 1.39 with parameter settings used in this 
paper), (b) A sensory blind zone oriented opposite the direction of forward motion of the particle 
with angular range 6. Particles are represented by small triangles. In this example, particles 6, 
8 and 9 are within particle I's blind zone, so their indices are not included in the set Si^t when 
generating repulsive/attractive forces acting on particle 1. 



sensory blind zones. Consequently, we revised the equations of motion using discrete time 
steps. The difference equations used for numerical simulation are 



m- 



At 



*>* _ /„, /3I„. |2 



At 

UMx) 



(5) 

(6) 
(7) 



where Xi^t and Vi,t are the position and the velocity of the i-th particle at time t, respectively; 
^i^t a randomly oriented vector with length 1 whose orientation is independent for each 
evaluation; Ui^t{x) the interaction potential surface for the i-th particle at time t; and Si^t 
the set of indices of all the particles whose positions are outside the blind zone of the z-th 
particle at time t (Fig. [D (b)). 
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We conducted numerical simulations of this model to produce a milling pattern similar 
m structure to those witnessed in schools of teleost fish, insects, and microorganisms [2], ^ 



0, y, 11, Q, Q- Specific values of fixed parameters are as follows: m = 1.0, Cr = 1.0, 
Ca = 0.5, Ir = 0.5, = 2.0, a = 1.6, (3 = 0.5. Initial conditions of each simulation were such 
that particles were randomly distributed within a square area of side length 2/^, and each 
particle was randornly oriented with magnitude of velocity randomly chosen from [0,feg] as 



described in 



111 



The model equations were numerically simulated from t = to 200 at interval At = 0.01. 
No spatial boundaries were enforced. For simulations recording the effect of stochastic force, 
7 was varied from to 10 at interval 0.5, while ^ = 0. For simulations testing the effect 
of sensory blind zones, 6 was varied from to 0.27r at interval O.Olvr, while 7 = 0. Each 
parameter setting was simulated using several population sizes, = 200, 300, 400, and 500. 
Ten simulation runs were conducted for each condition. 

Several metrics were used to characterize the simulation results. These include average 
absolute velocity Vabs, ratio of halting particles H, normalized angular momentum M, and 
normalized absolute angular momentum Mabs, defined as follows (same or similar metrics 
were used in jol. [l^): 

Vabs = El^^l/^ (8) 

i 

H = S.t. \Vi\ <flVe,}\/N (9) 

M = (10) 



{h 


s.t. 

n > 


< Vi 
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Vi\ 






Vi\ 



Mabs = t:'; ,, ' (11) 



Here = — Xc where Xc is the swarm's center of mass. To measure H we used 20% of the 
equilibrium velocity (/i = 0.2) as a threshold to determine whether a particle was halting 
or not. When used comparatively, M and Mabs make it possible to distinguish single mill 
formation from double mill formation in which two mills rotate with opposite sense around 
similar, but not identical, centers of mass j^, [ll|. All the metrics were averaged over the 
last 10 time steps of each simulation. 

Figs. [2] and [3] depict the processes of structural decay produced by stochastic force and 
blind zone perturbations on the milling behavior of 500 particles. A transition from milling 
state to disordered state was induced by increasing the magnitude of stochastic force across 
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7 ~ 7.0 (Fig. [2] (a), Fig. [3] (a)-(c)). When the range of sensory bhnd zones 6 was increased, 
however, structural degradation was very different, involving a new spatially heterogeneous 
transition (Fig. [2] (b). Fig. [3] (d)-(f)). Initiation of collapse occurred at ^ O.OStt, where 
particles near the center of the mill ceased rotation and formed a stationary core that has not 
previously been described. We call this new state a "carousel" state. This state is different 
from the rigid-body rotation reported in [o], [lli and the compact but disordered state of [s] 
because it is characterized by a sharp boundary between the milling surface and the central 
stationary core made of particles with near zero velocity (Fig. [3] (e)). A secondary transition 
was observed across 9 ~ O.OGvr where the particles moving in the periphery became abruptly 
disordered and lost coherence in motion while the particles in the central core remained 
stationary (Fig. [3] (f)). We call this concluding state a "surface disordered" state. 

Particles inside the core in the carousel and surface disordered states lose their velocity 
due to a perceptual and consequent force asymmetry. A sensory blind zone creates a lon- 
gitudinal imbalance between forces derived from particles ahead and forces from particles 
behind. Because the pairwise equilibrium distance Teq is much longer than the character- 
istic distance between neighboring particles in a swarm, the imbalance takes effect in the 
regime of repulsive interactions and thus results in a net resistance against self-propulsion 
of particles. A particle near the center of the swarm has more particles to its front and 
back than does a particle rotating in the periphery, since the density of particles is inversely 
related to the distance from the center of the swarm under the parameter set used here 10 1 
(numerically confirmed in our simulation results; data not shown). Thus, the net resistance 
against forward motion resulting from the blind zone is larger for particles rotating close to 
the mill's center. Within a certain distance to the mill's center, the resistance exceeds the 
range of self-propulsive force possible in Eq. ([2]), and consequently particles cease motion. 
This transition does not occur due to increasing 7 because the effect of stochastic force is 
spatially isotropic: it is equally likely to force a particle in any direction and therefore does 
not lead to cessation of movement. 

Fig. m summarizes all the simulation results, showing the dependence of the final values of 
Vabsi H, M and Mabs on 7, 6 and A^. The onset and the mechanism of structural degradation 
of milling behavior are different between cases with increasing 7 and 6. The milling structure 
is fairly robust to small 7, and it suddenly collapses at 7 ~ 7.0, nearly independently of A^. 
The H plot shows no particles halting in this transition. In contrast, plots of increasing 9 
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FIG. 2: A visual comparison of the effects of increasing stochastic force and the effects of increasing 
range of sensory bUnd zones on the milhng behavior of 500 particles. Each image is a final snapshot 
of a simulated particle swarm taken at t = 200. Particles have tails that represent the orientation 
and magnitude of their velocity, (a) Results with increasing stochastic force 7 while 9 = 0.0. 
Transition from milling to disordered states occurred at 7 ~ 7.0. (b) Results with increasing range 
of sensory blind zones 6 while 7 = 0.0. Transitions from milling to carousel and from carousel to 
surface disordered states occurred at f» O.OSvr and 9 Rr: O.OGtt, respectively. See also Fig. [3l 
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FIG. 3: Tangential velocities of particles at a distance r from the center of mass. Data were 
obtained from numerical simulations of 500 particles at t = 200. The direction of rotation of the 
majority was taken as positive. 

illustrate that structural degradation in increasing is a two-fold process. The first transition 
from milling to carousel was detected in Vabs and H (emergence of halting particles and 
consequent decrease of average velocity). The second transition from carousel to surface 
disordered was detected in M and Mats (loss of coherence in angular momentum). It was 
also observed in our results that the onsets of these transitions depended significantly on 
A^. This can be understood in that larger increased the density at the mill's center and 
made particles more reactive to blind zone induced halting. 

The blind zone ranges used in these simulations were extremely small from a biological 
viewpoint. The largest value tested, 9 = 0.2n, was just 10% of the perception range, which 
was more than sufficient to destroy milling patterns in all cases. Kinematic SPP models, 
on the other hand, can produce and maintain milling patterns despite considerable sensory 
blind zones jl4]. To understand this discrepancy in model properties, we note one important 
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FIG. 4: (Color online) Comparison of the values of the four metrics (average absolute velocity Vabs, 
ratio of halting particles H, normalized angular momentum M and normalized absolute angular 
momentum M^bs) measured for over all simulations for N = 200, 300, 400, 500. Each data point 
represents the average of ten simulation runs with an error bar, measured in standard deviations, 
(a) Results with increasing stochastic force, where the collapse of milling behavior is always reached 
at 7 7.0 regardless of N. (b) Results with increasing range of sensory blind zones. The collapse 
of mill behavior is twofold; the top two plots {Vabs and H) capture the first transition from milling 
to carousel, while the bottom two (M and Mabs) capture the second transition from carousel to 
surface disordered. The critical values of 9 for these transitions depend on N, indicating that larger 
populations are more susceptible to blind zone perturbations. 

9 



difference between Icinematic and dynamic frameworks: while particles in kinematic models 
are always constrained to move with a non-zero velocity, there is a possibility for particles 
to halt in dynamic models that becomes significant in the presence of rear blind zones like 
those assumed in our model. This leads us to a hypothesis that milling behavior in an ag- 
gregate of organisms may sensitively depend on their ability to maintain constant velocity. 
Specifically, for organisms that keep moving autonomously at a near constant pace, milling 
behavior emerges relatively easily even with considerable sensory blind zones. In contrast, 
for organisms whose motion strongly depends on (either sensory or physical) environmental 
stimuli, milling behavior requires a nearly complete panoramic range of interaction, espe- 
cially to perceive the pressure from behind and gain enough forward propulsion to maintain 
constant velocity. 

There are several biological observations that directly or indirectly support our hy- 
jothesis. Milling behavior is often reported in groups of microorganisms and insects 
2I, 3, y, 0, Is, ^, 10, 19]- These organisms rely chiefly on direct physical contact and chemical 
sensory input, respectively, when forming mills and therefore use omnidirectional sensory 
capabilities to form aggregates. They also have the ability to cease motion. Additionally, 
there is a great deal of evidence supporting the isotropic sensory ability of the lateral line 
in teleost fish, some of which demonstrate milling behavior |20|. A recent study that in- 
vestigated the superficial organization of neuromasts composing the lateral line in goldfish 
showed that neuromasts' most sensitive axes were oriented in almost every direction 2ll |. 
Moreover, a recent study on Mormon crickets 22[ reports that the physical, cannibalistic 
threat of protein and salt deprived individuals from behind plays a critical role in creating 
a large-scale coherent march. When some crickets are immobilized and therefore unable to 
respond to a push from behind, the march halts. This study provides clear evidence support- 
ing our conjecture that inputs (physical pressure) from behind a particle can be important 
in the formation of coherent swarming patterns. 

In summary, we computationally studied the effects of sensory blind zones on the stability 
of self-organizing mill formation in a dynamic SPP model. We found that milling behavior 
collapses through two spatially distinct transitions in response to an increasing range of rear 
blind zones, characterized by a halting regime emanating from the center of the swarm and 
then a disorganization of coherent motion in the periphery area. This is quite different from 
pattern collapse observed with increasing stochastic force described by a spatially uniform 
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transition to a compact but disordered state [5|]. Combined with other results obtained with 
kinematic SPP models, our results suggest a possible relationship between collective behavior 
and sensory systems of biological organisms: species that engage in mill formation in nature 
may necessarily have an omnidirectional sensory system if they do not maintain constant 
velocity by themselves. This is a hypothesis testable and falsifiable through experimental 
observation. 
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